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ABSTRACT 

We present images from a space- VLBI survey using the facilities of the VLBI Space Observatory 
Programme (VSOP), drawing our sample from the well-studied Pearson- Readhead survey of extragalactic 
radio sources. Our survey has taken advantage of long space- VLBI baselines and large arrays of ground 
antennas, such as the Very Long Baseline Array and European VLBI Network, to obtain high resolution 
images of 27 active galactic nuclei, and to measure the core brightness temperatures of these sources 
more accurately than is possible from the ground. A detailed analysis of the source properties is given 
in accompanying papers. We have also performed an extensive series of simulations to investigate the 
errors in VSOP images caused by the relatively large holes in the (it, v) plane when sources are observed 
near the orbit normal direction. We find that while the nominal dynamic range (defined as the ratio 
of map peak to off-source error) often exceeds 1000:1, the true dynamic range (map peak to on-source 
error) is only about 30:1 for relatively complex core-jet sources. For sources dominated by a strong point 
source, this value rises to approximately 100:1. We find the true dynamic range to be a relatively weak 
function of the difference in position angle (PA) between the jet PA and (u, v) coverage major axis PA. 
For low signal-to-noise regions typically located down the jet away from the core, large errors can occur, 
causing spurious features in VSOP images that should be interpreted with caution. 

Subject headings: galaxies : jets — galaxies : active — quasars : general — radio galaxies : continuum 
— techniques: high angular resolution — techniques : interferometric 



elliptical orbit has an apogee height of 21,400 km, mak- 
ing it possible to directly image compact radio sources at 
1.4 and 5 GHz with sub-milliarcsecond angular resolution, 
and to measure brightness temperatures in excess of 10 12 
K (Hirabayashi et al. 1998). 

The images presented here are of the highest angular 
resolution ever achieved at 5 GHz, and reveal substan- 
tially more structural detail than previous ground VLBI 
studies of the Pearson- Readhead sample (e.g., Pearson & 
Readhead 1981, 1988; Cawthorne et al. 1993). In §2 we de- 
scribe our VSOP observations and data analysis methods, 
and present the images and basic data on our sample. We 
describe our model fitting technique in §2.2, and discuss 
the impact of incomplete (u, v) plane coverage on space- 
VLBI images in §3. A full analysis of our space- VLBI data 
is presented in two companion papers (Tingay et al. 2001; 
Lister et al. 2001). 

2. OBSERVATIONS AND DATA ANALYSIS 

The Pearson-Readhead (PR) survey is a complete sam- 
ple of 65 northern (S > +35°) AGNs with 5 GHz flux 
densities > 1.3 Jy, and galactic latitude \b\ > 10° (Pear- 
son & Readhead 1988). In order to maximize the likeli- 
hood of detecting fringes to HALCA (at the 7a level), we 
have selected as space- VLBI targets 31 PR objects that 
have flux densities > 0.4 Jy on the longest Earth baselines 
(~ 10,000 km). 

Observations of 27 of these objects were conducted with 

l Jet Propulsion Laboratory, California Institute of Technology, MS 238-332, 4800 Oak Grove Drive, Pasadena, CA 91109-8099 
2 Current address: NRAO, 520 Edgemont Road, Charlottesville, VA 22903-2454 
3 Australia Telescope National Facility, P.O. Box 76, Epping, NSW 2121, Australia 
4 Whittier College, Department of Physics, 13406 Philadelphia Street, Whittier, CA 90608 
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1. INTRODUCTION 

Ground-based very long baseline interferometric (VLBI) 
surveys of the parsec-scale regions of active galactic nuclei 
(AGNs) have led to important insights into the nature of 
relativistic bulk outflows in compact radio sources. Here 
we present observational data from the first imaging sur- 
vey of a complete sample of compact extragalactic radio 
sources made with space- VLBI. Our objects are selected 
from the well-studied Pearson-Readhead sample of AGNs, 
and are representative of some of the brightest, most com- 
pact AGNs in the sky. 

Space- VLBI represents a logical step in the evolution of 
radio interferometry, and has the potential to provide the 
highest angular resolution images of any observing tech- 
nique. Since our understanding of AGNs has historically 
been closely tied to resolution (e.g., the discovery of jets 
with connected interferometers, and relativistic bulk mo- 
tions with VLBI), space- VLBI has an important role in 
AGN research. 

The first successful space- VLBI experiments were ac- 
complished in 1986 and 1987 using a TDRSS communi- 
cations satellite (Levy et al. 1986; Linfield et al. 1989, 
1990), and provided direct evidence for exceedingly com- 
pact emitting regions in extragalactic radio sources. In 
1997, the first satellite dedicated to radio astronomical 
observations (HALCA) was launched by Japan's Institute 
of Space and Astronautical Science, as part of the VLBI 
Space Observatory Programme (VSOP). HALCA's highly 



2 



Lister et al. 



the HALCA orbiting antenna over the period 1997 Au- 
gust to 1999 April at a frequency of 5 GHz, in conjunc- 
tion with telescopes from either the Very Long Baseline 
Array (VLB A), the European VLBI Network (EVN) , or 
in some cases, both arrays. Four sources (0454+844, 
0804+499, 0945+408, and 0954+658) were not observed 
due to scheduling constraints. HALCA was tracked typi- 
cally between four and five hours per observation, synthe- 
sizing (u,v) coverages with a maximum baseline of nearly 
500 MA, or 30, 000 km. A journal of our observations is 
given in Table 1. 

Left circular polarization data from the ground anten- 
nas, or in the case of the orbiter, the spacecraft track- 
ing station, were recorded on VLBA format tapes in two 
16 MHz IFs and shipped to the VLBA correlator in So- 
corro for processing. The correlated data were amplitude- 
calibrated, fringe-fitted, and frequency-averaged in AIPS 5 
before being exported to DIFMAP (Shepherd et al. 1994) 
for self-calibration, imaging and model-fitting. 

Once imported into DIFMAP, the visibilities were edited 
and then coherently averaged over 30 second intervals. 
Due to its relatively small diameter (8 m) and poor an- 
tenna efficiency, the system equivalent flux density of 
HALCA is approximately 50 times worse than that of a 
typical VLBA antenna at 5 GHz. In order to give effec- 
tively equal weight to the ground-only and ground-space 
visibilities, we increased the weight of the HALCA an- 
tenna in self-calibration and subsequent model-fitting by 
a factor based on the relative number of ground-ground 
and space-ground visibilities and their respective baseline 
sensitivities. These factors are listed in Table 1. 

Following amplitude and phase self-calibration, we pro- 
duced several images for each source using different visi- 
bility weighting schemes. The choice of weighting scheme 
has a large influence on space- VLBI images, due to the 
unequal numbers of ground-ground and ground-space vis- 
ibilities, and the large differences in baseline noise level 
described above. Uniform weighting reduces the influ- 
ence of the ground-ground visibilities on the image, which 
would normally dominate since they outnumber the space- 
ground visibilities. In the right-hand panels of Figures 1-9, 
we show uniformly-weighted images made using all avail- 
able baselines. These images have no additional weighting 
by amplitude errors. In the left-hand panels, we show 
naturally-weighted images made using ground-ground vis- 
ibilities only. We also include a representation of the 
restoring beam in the lower left corner of each image, and 
a bar showing the spatial scale. The latter is based on 
a standard Frcidmann cosmology with zero cosmological 
constant, deceleration parameter q a = 0.1, and Hubble 
constant H Q = 100 h km s^ 1 Mpc -1 . We will assume 
this model throughout this paper. The parameters of the 
restoring beam, the peak brightness, and the contour levels 
of each image are listed in Table 2. 

In Table 3 we list the basic data on each source, in- 
cluding redshift, total amount of cleaned flux in the VLBI 
map, and an estimate of the total (single-dish) 4.8 GHz 
flux density at the epoch of the VSOP observation. The 
latter data were obtained by interpolating data from the 

5 ALPS and the VLBA are maintained by the National Radio As 
under cooperative agreement with the National Science Foundation. 

6 http://www.astro. lsa.umich.edu/obs/radiotel/umrao. html 



University of Michigan monitoring program 6 , and have an 
approximate error of ±0.1 Jy. 

2.1. Visibility Functions 

The visibility function of a source represents the ob- 
served power measured on different spatial scales, with 
the long baselines representing small-scale structure, and 
vice- versa. In Figure 10 we plot the upper envelopes of the 
visibility amplitude distribution versus projected baseline 
length for the strongest and weakest sources, respectively. 
Instead of dropping smoothly to zero, the curves show a 
pronounced flattening at baselines greater than an earth 
diameter. This is indicative of strong, compact compo- 
nents that remain unresolved at the longest baselines. The 
brightness temperatures of these components and their 
implications for relativistic beaming models are discussed 
more fully by Tingay et al. (2001). 

2.2. Component Model- fitting 

We began our model-fitting analysis of each source with 
the core component, which is commonly believed to be as- 
sociated with a flat-spectrum component located near the 
base of the jet. Wherever possible, we relied on previously 
published identifications to determine the position of the 
core. For sixteen sources the core (labeled "D" in Fig- 
ures 1-9) was located at the end of the jet, and was also 
the brightest component in the 5 GHz image. In the case 
of nine sources, the component at the end of the jet was 
the brightest component in images at higher frequencies, 
but not at 5 GHz. The core component identifications for 
the two radio galaxies in our sample are somewhat un- 
certain, as their jet structure may not be one-sided. For 
3C 84 (0316+413), we arbitrarily assigned the core posi- 
tion to the northernmost component of the jet, while for 
2021+614, we used the core identification of Tschager et 
al. (2000). 

We first removed the set of point source clean compo- 
nents used to fit the core from the clean component model 
and replaced them with an elliptical Gaussian described by 
6 free parameters. We then fit the Gaussian in the (u, v) 
plane using the MODELFIT task in DIFMAP, which uses 
a linear least squares fit to the amplitude and phase to find 
the best fit model. For all our fits, we adjusted the visibil- 
ity weights used in the least squares calculation, such that 
the ground antennas and the space antenna contributed 
equally to the reduced chi-squared statistic. We inter- 
spersed runs of MODELFIT with phase self-calibration, to 
ensure that the reduced chi-squared was minimized to find 
the best fit models. The parameters of our core-component 
fits are given in Table 4. 

After fitting the core, we removed the clean components 
associated with various other isolated features and model- 
fit them with elliptical Gaussians. In order to avoid in- 
troducing new nomenclature for individual components, 
we have used existing names from the literature wherever 
possible. These are given in Table 5, along with the pa- 
rameters of our fits. We caution that these fits may not be 
robust, especially in regions where there is continuous jet 
emission. Indeed, some features (e.g., the southernmost 

ny Observatory, which is operated by Associated Universities, Inc., 
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component of 2021+614) were too complex to be repre- 
sented by simple Gaussians, and were not model-fitted. 

We determined errors on the parameters of the core- 
component models using the DIFWRAP package (Lovcll 
1998), a graphical interface to DIFMAP which fixes given 
parameters of the models at different values which are per- 
turbations on the best-fit values, and allows all other pa- 
rameters to vary freely to re-minimize the reduced chi- 
squared. In this way it is possible to determine how far 
from the best-fit value a given parameter can be forced be- 
fore the model no longer fits the (u, v) plane data, defining 
the error on that parameter. 

For each source, we first determined whether the core 
component could be fit with an unresolved point source 
(i.e., a delta- function component). If a point source did 
not fit the visibilities within the amplitude errors, we then 
tried one-dimensional (zero-axial ratio) components of var- 
ious lengths and position angles (PAs). If none of these 
provided a good fit, we varied the parameters of the best- 
fit model as described above to determine the range of 
possible brightness temperature that still fit the visibility 
data. In two cases (0923+392 and 2021+614), the cores 
were too weak with respect to the extended emission to 
determine a best fit range of brightness temperature. 

The components in our sample could have also been fit 
using the surface brightness distributions of optically thin 
or optically thick spheroids, instead of Gaussian profiles. 
This is due to the fact that their Fourier transforms are 
virtually identical in the (u,v) plane out to the point at 
which the flux density is 50% of the zero baseline flux 
density (Pearson 1995). Since we do not generally observe 
past this point for the core components in our sources, we 
have adopted Gaussian surface brightness distributions. 
To convert the observed brightness temperature assuming 
a Gaussian component to an optically thin or optically 
thick sphere, corrective factors of 0.67 and 0.56 respec- 
tively, should be applied. 

We find that the sizes of our fitted components gen- 
erally increase exponentially with distance down the jet. 
In Figure 11 we present a log-log plot of fitted compo- 
nent area (ft) in square parsecs against projected distance 
from the core component (r) in parsecs. A strong trend is 
present, with a linear correlation coefficient of 0.42, and a 
significance level of 99.985%. The dashed line represents 
a linear regression through the points (in log space) of the 
form logft = (1.4 + 0.1) logr+ (0.07 + 0.1). The fit consid- 
ering quasars alone is logSl = (1.1 + 0.2) logr + (0.5 + 0.2), 
and for BL Lacertae objects alone it is logfi = (1.5 ± 
0.2) logr - (0.2 ± 0.2). These fitted slopes are all signif- 
icantly less than two, which is the value for a conical jet 
of constant opening angle. Considering the fact that the 
projected distances in Fig. 11 are likely to be severely fore- 
shortened due to small viewing angles (see Paper II), the 
trend indicates that the jets are still in a state of expan- 
sion and are undergoing collimation at distances of several 
hundred parsecs from the core. 

2.3. Influence of (u,v) Coverage on Model-fits 

The (u, v) coverages of most of our space- VLBI obser- 
vations are highly elongated, due to the elliptical orbit of 
the HALCA satellite (Figure 12). As a result, the angular 
extent of the jets and components in our sample are much 



better constrained along one direction. We have examined 
the possible effects of this unequal (u, v) plane sampling on 
our core component model-fits by inter-comparing the po- 
sition angles of the fitted Gaussian component major axis, 
the position angle of the restoring beam, and the position 
angle of the parsec-scale jet. 

We find that the major axis position angle of the fit- 
ted core component is strongly correlated with that of the 
beam and (u, v) coverage, and therefore does not likely re- 
flect an intrinsic property of the source. For those sources 
whose cores could be fit with a zero-axial ratio (ZAR) com- 
ponent, the (u, v) coverage is generally perpendicular to 
the jet direction, and the fitted major axis of the ZAR com- 
ponent is aligned with the jet. This suggests that there are 
two or more closely spaced components in the core region 
that are separated by much less than a beam width. Since 
they cannot be resolved, the visibility data can be fit with 
a long, narrow (ZAR) component. In this case, the upper 
limit on the brightness temperature is unconstrained. 

3. EFFECTS OF INCOMPLETE (U, V) PLANE COVERAGE 
ON SPACE- VLBI IMAGES 

Most of the space- VLBI observations of our sample were 
scheduled when the source lay close to the orbit normal. 
Under these conditions the maximum possible spatial res- 
olution can be achieved, at the cost of creating relatively 
large holes in the (u, v) plane aperture coverage. This 
is due to the the fact that the HALCA spacecraft has 
a highly elliptical orbit with an apogee height of 21,400 
km, which greatly exceeds an Earth diameter. The ef- 
fects of (u, v) plane under-sampling are well-documented 
for ground-based inteferometric arrays (e.g., Perley 1989), 
but relatively little work has been published for the case 
of space- VLBI. Here we discuss the effects of (u, w)-holes 
on the quality of our space- VLBI images, and how the re- 
sulting imaging errors depend on the source's jet position 
angle relative the position angle of the major axis of the 
(u, v) coverage. 

3.1. Example Simulations 

We have addressed these issues by performing exten- 
sive image simulations using several different source mod- 
els. These simulations have thermal amplitude and phases 
errors added to them at the same level as is present in 
the real data. Thus, the RMS noise-level on simulated 
HALCA baselines is about 7 times higher than on the 
simulated ground-only baselines. The simulated data were 
imaged using the same DIFMAP package that was used to 
reduce our actual space- VLBI data. For the purposes of 
this discussion we will consider the image artifacts which 
are present in 5 GHz images for the following source mod- 
els, all of which are composed of Gaussian components: 

1) A complex core-jet source which has a total flux den- 
sity of 2.7 Jy. We will refer to this model as cj-complex; 

2) A simple core-jet source, which has a 1 Jy, 10 12 K 
core with axial ratio unity and a 1 Jy, 10 11 K 'jet' with 
axial ratio 1/3 whose centroid is located at a distance of 
1.21 milliarcseconds away from the core. The major axis 
of this component measures 1.21 milliarcseconds, and is 
aligned with the core-jet direction. We will refer to this 
model as cj-simple; 

3) A 1 Jy, 10 12 K circular Gaussian; 
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4) A 1 Jy, 10 14 K circular Gaussian. This source is com- 
pletely unresolved for VSOP 5 GHz observations and can 
be used to determined the image artifacts in unresolved 
point sources. It is therefore a useful point of reference 
when examining the artifacts in more complex sources. 

In Figures 13 and 14 we show two simulations we have 
performed using the cj-complex and cj-simple models re- 
spectively. We use the same (u, v) coverage in both cases 
(panel b), which has a major axis PA of 0°. Both jets have 
a PA of 60°. 

Panels (a), (c), and (d) show the model, simulated im- 
age, and difference image, respectively, all convolved with 
an identical beam that corresponds to uniform weighting 
of the visibility data. The difference image is simply de- 
fined as the difference between the the simulated image 
and the model. The difference images reveal a common 
error feature seen in the imaging of core-jet sources, in the 
form of a sinusoidal ripple that runs along the length of the 
jet. This feature is the result of a well-known instability 
of the CLEAN algorithm, and has modulations that tend 
to correspond with the spacing of holes in the (u, v) plane 
(Cornwell 1983). In our simulations, this error feature is 
so prominent that it can cause visible artifacts to occur 
in the jet that may easily be over-interpreted as discrete 
components. It is important, therefore, to exercise caution 
in the interpretation of relatively weak features in VSOP 
images. 

A good example of this can be found in the model in Fig- 
ure 13, in which panel (a) shows a rather smooth jet. The 
simulated VSOP image, on the other hand, shows a series 
of discrete components which run down the jet (panel c). 
In panels (d) and (f ) the difference map signal-to-noise ra- 
tio (SNR) at each pixel and the percentage error in the 
simulated are shown as a function of the simulated image 
SNR. The SNRs quoted here use a reference noise level de- 
termined far away from the source structure. From these 
two panels we can see how the imaging errors decrease as 
the map SNR increases. For high SNR points (>100) the 
image error is typically less than 10%. However, for in- 
termediate SNR points (10 < SNR < 100) the percentage 
error can almost be 100%, so caution must be exercised in 
over-interpreting data with this SNR. 

3.2. Dynamic Range 

A common method of quantifying image errors is to take 
the maximum in the image divided by the off-source RMS 
noise level (DR^ms), which is sometimes referred to as the 
dynamic range of an image. For example, the simulated 
images shown in Figures 13 and 14 have D-Rrms values 
of 1600:1 and 2800:1 respectively. However, as our simula- 
tions show, this is not the true dynamic range in the image 
as the on-source errors are much larger than the off-source 
errors. We define here a new dynamic range, called -D-Rdifi 
which is the maximum in the image divided by the maxi- 
mum value (whether positive or negative) in the difference 
image. Using this definition, DRdie has values of 34:1 and 
28:1 for the images in Figures 13 and 14. Consequently the 
true dynamic ranges (DR^ie ) of these images are a factor 
of ~ 50 worse compared to the dynamic range using the 
off-source noise level (D-Rrms)- For comparison, the im- 
ages made using the single 10 12 K and 10 14 K Gaussian 
component models, have values of DRrms of 3000:1 and 



4000:1, and DR dm of 48:1 and 130:0, respectively. Thus, 
we see a general trend of increasing image fidelity as the 
source structure becomes simpler (i.e., DR^m decreases). 
However, even in the limit of a unresolved point source, 
the value of the true dynamic range is limited to approxi- 
mately 100 : 1. 

3.3. Effects of Jet Position Angle on Image Dynamic 

Range 

One might expect that the dynamic range that we can 
obtain on a given observation of a source should depend 
on the difference between a source's jet PA and the PA of 
the (u,v) coverage major axis (APA). For APA =0° we 
get the highest resolution along the jet and for APA = 
90° we get the highest resolution perpendicular to the jet. 
Thus, the ability of the sinusoidal image error to develop 
could depend on the value of APA. Consequently, certain 
alignments of jet position angle with respect to the (u, v) 
coverage might be more reliable than those of a different 
APA. We have investigated this phenomenon by conduct- 
ing a series of simulations using the same (u, v) coverage 
as shown in Figures 13 and 14 but for different values of 
the jet PA in the model sources cj-complex and cj-simple. 
In Figures 15 and 16 we plot how both D-Rrms and -D-Rdiff 
vary as a function of the jet PA while keeping the PA of 
the major axis of the (u, v) coverage fixed. As can be seen, 
DRrms is almost independent of APA and DR^is is only 
a relatively weak function of this position angle difference. 

4. SUMMARY 

We have presented data from the first imaging survey of 
a complete sample of compact extragalactic radio sources 
made with space- VLBI. These data represent the highest 
angular resolution images ever obtained for these objects 
at 5 GHz, and provide important information regarding 
the structure of AGN jets on parsec scales. In particular, 
we have found an exponential trend of increasing compo- 
nent size with distance down the jet, which indicates that 
these jets are still undergoing expansion on scales of sev- 
eral hundred parsecs. We have also detected a pronounced 
flattening in the visibility functions of our sources, which 
is indicative of strong, compact components that remain 
unresolved at the longest space-Earth baselines. We dis- 
cuss the implications of this trend and the other general 
properties of our sample in two companion papers (Tingay 
ct al. 2001; Lister ct al. 2001). 

We have performed an extensive series of simulations to 
investigate the image errors in VSOP images caused by 
the relatively large holes in the (u, v) plane when sources 
are observed near the orbit normal direction. We find that 
while the nominal dynamic range (-D-Rrms) often exceeds 
1000:1, the true dynamic range (Di?diff ) is about 30:1 for 
relatively complex core-jet sources. For sources dominated 
by a strong point source, this value rises to approximately 
100:1. The true dynamic range is also found to be a rel- 
atively weak function of the difference in position angle 
between the jet PA and (u, v) coverage major axis PA. 
For high SNR regions in the image (SNR>100) the error 
on individual pixels is typically less than 10%. However, 
for low SNR regions typically located down the jet away 
from the core, large errors can occur and spurious features 
can be seen. Caution should therefore be exercised when 
interpreting regions of low SNR in VSOP images. 
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Table 1 
Journal of Observations 



IAU 


Other 


Observing 


Ground 


Observing 


Integration 


HALCA 


Name 


Name 


Date 


Antenna Array 


Freq.[GHz] 


Time [h] a 


Weight b 


0016+731 




1998 Mar 2 


VLBA+U 


4.81 


10 (3) 


1156 


0133+476 


OC 457 


1999 Aug 16 


VLBA 


4.81 


5 (4.5) 


274 


0153+744 




1998 Aug 18 


VLBA-Br,Sc,La 


4.81 


6 (5.5) 


187 


0212+735 




1997 Sep 5 


VLBA-Sc 


4.97 


8 (4.5) 


415 


0316+413 


3C 84 


1998 Aug 25 


VLBA+Y 


4.97 


6 (5) 


538 


07H+356 


01 318 


1999 Apr 9 


VLBA 


4.81 


r> / a \ 

6 (4) 


324 


0814+425 


OJ 425 


1999 Apr 24 


VLBA 


4.81 


6 (5) 


289 


0836+710 


4C 71.07 


1997 Oct 7 


VLBA 


4.97 


5.5 (3.5) 


390 


0859+470 


4C 47.29 


1999 Feb 14 


JEGNMTWO 


4.97 


6.5 (1) 


1460 


0906+430 


3C 216 


1999 Feb 14 


JEGNMTW 


4.97 


5.5 (3) 


399 


0923+392 


4C 39.25 


1997 Oct 23 


VLBA+Y 


4.97 


2.5 (2) 


847 


1624+416 


4C 41.32 


1998 Feb 7 


VLBA-Mk 


4.81 


4 (1.5) 


347 


1633+382 


4C 38.41 


1998 Aug 4 


VLBA+Y 


4.81 


5.5 (4) 


531 


1637+574 


OS 562 


1998 Apr 21 


VLBA-Sc,Ov,Kp 


4.81 


6.5 (4.5) 


280 


1641+399 


3C 345 


1998 Jul 28 


VLBA+EY 


4.81 


4.5 (3.5) 


674 


1642+690 


4C 69.21 


1998 May 31 


JEGNMWS 


4.97 


6.5 (4.5) 


225 


1652+398 


MK 501 


1998 Apr 7 


VLBA+E 


4.81 


5 (3.5) 


370 


1739+522 


4C 51.37 


1998 Jun 14 


VLBA 


4.81 


6.5 (4) 


391 


1749+701 




1999 Jun 1 


VLBA 


4.81 


8(5) 


406 


1803+784 




1997 Oct 17 


VLBA 


4.97 


6(4) 


370 


1807+698 


3C 371 


1998 Mar 11 


VLBA 


4.81 


6 (5) 


349 


1823+568 


4C 56.27 


1998 May 31 


JEGNMO 


4.97 


4.5 (3.5) 


542 


1828+487 


3C 380 


1998 Jul 4 


VLBA+E 


4.81 


6(4) 


584 


1928+738 


4C 73.18 


1997 Aug 22 


VLBA+E 


4.97 


5.5 (4.5) 


482 


1954+513 


OV 591 


1997 Nov 10 


EGNMWO 


4.97 


12 (4) 


777 


2021+614 


OW 637 


1997 Nov 7 


VLBA+ENMO 


4.97 


5.5 (4.5) 


575 


2200+420 


BL Lac 


1997 Dec 8 


VLBA 


4.97 


6.5 (5) 


315 



Note. — Antenna abbreviations: Br = Brewster, E = Effclsberg, G = Green Bank, J=Jodrell 
Bank, Kp = Kitt Peak, La = Los Alamos, M = Medicina, N = Noto, O = Onsala, Ov = Owens 
Valley, S = Sheshan, Sc = Saint Croix, T = Torun, U = Usuda, W = Westerbork, Y = Phased 
Very Large Array. 

a Approximate ground integration time on source. Values in parentheses indicate integration time 
that included the HALCA spacecraft. 

b Wcighting assigned to the HALCA antenna during self-calibration and model-fitting (ground 
visibilities had a nominal weighting factor of unity). 
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Table 2 

Summary of Image parameters 



Source 
(1) 


Baselines 
(2) 




Beam 
(3) 




PA 

(4) 


Peak 
(5) 


Contour levels 
(6) 


nm ft i 7qi 
UU1D-|- (ol 


Ground 


1 
1 


1 


X 


9 

z 


H9 
Uz 


— lo 


n 
U 


zi^n 

40U 


n nfi 
— u.uo, 


£1 44 
. . . , D1.44 




A 11 

All 


U 


38 


X 





77 
< f 


■7ft 
— ( O 





216 


1.0, . 


96 


n 1 QQ_L /I 7ft 


Ground 


1 


79 
/ z 


X 




HQ 

uy 


Q 
O 




f yo 


n 1 

— u.i, . 


51 2 




A 11 

All 


U 


23 


X 





01 


OO 

— zy 


I 


340 


—0.5 


64 


Uloo-f- ( 44 


Ground 


1 


uo 


X 


9 

z 


1 n 

1U 


on 
— OU 


n 
U 


9Q^ 

zyo 


1 

— 1, . . . 


, 04 




A 11 

All 


n 
U 


36 


x 





OZ 


Qft 
— OO 


Q 


070 


—A 
4, . . . 


64 

, D4 


HOI O I 7Q c 

Uzlz4- ( 00 


Ground 


1 
1 


fi^ 
00 


X 


9 

z 


O 1 

y4 


7ft 
( O 


9 
Z 


^^7 
oO 1 


n 1 
— u.i, . 


^1 9 
. . , 01. z 




A 11 

All 


n 
U 


24 







ft7 

f 


Oo 





786 


— n 4 

U.4, . 


51.2 


n q i ft _i_ /i 1 q 


Ground 


1 
1 


97 

Z I 


x 


2 


sn 

oU 


9ft 

— zo 




944 


— 1.0, . 


. . , yo 




A 11 

All 


U 


31 


x 





77 
< / 


n 
U 





375 


— 9 ^ 
z.o, . 


80 


U / 11 + OOD 


Ground 


1 
1 


OO 


X 






Oft 

zo 


A 

— 4 


n 
u 


7n7 

1 u 1 


n 9 


M 9 

. . , 01. z 




A 11 

All 


n 
U 


01 


x 


n 
u 


00 


ftQ 
-DO 


n 

u 


189 


—4 


, O^ 


nsi /i 1 /Iok 


Ground 


1 
1 


cc 



x 


■3 




"I I 


91 

— Zl 


n 
u 


ooy 


n 9 

— u.z, . 


51 2 




A 11 

All 


n 
U 


26 


x 





on 
oZ 


A O 

— 4z 





437 


—0.75 


96 


nooc 1 71 n 

U00O + ( 1U 


Ground 


1 
1 


^7 
/ 


X 


9 

z 


go 
OO 


K "I 

ol 


u 


044 


n 9^ 
— u.zo, 


(\A 
. . . , 04 




A 11 

All 


n 

u 


22 


x 


n 

u 


ftO 
OZ 


Qft 
— OO 


u 


^79 
/ z 


— 1 25 


p.n 

. . . , OU 


OQf;0_L/l "7fl 

uooy+4 ( u 


Ground 


A 

u 


.JO 


x 


2 


on 
zU 


1 O 

— iy 


u 


000 


n ^ 


fi4 




A 11 

All 


("1 

u 


20 




1 


og 

OO 


OO 

— OZ 





264 


—5 


80 


nnriR 1 /ion 

uyuo+4ou 


Ground 




z 


9"< 
zo 


x 


O 


TCI 


on 
— oU 


n 

u 


fii n 

O-LU 


n fi^ 


. . . , oo.z 




A 11 

All 


n 
U 


"1 7 


x 


n 
u 


on 
OZ 


A A 
— 44 


n 

u 


OOI 


— -L.O, . 


. . , yo 


UyzO-f-oyz 


Ground 


1 
1 


i 


x 




ftQ 
Do 


1 ft 
— 10 


9 


Ann 


n 1 


51 2 




A 11 

All 


n 
U 


98 
zo 


x 


2 


no 
zo 


1 Q 

— lo 




nQ9 
uyz 


n ^ 

— u.o, . 


R4 


1 ftO/1 1 /11ft 
I0z4-t-4l0 


Ground 


1 


on 


X 


4 

4 


-1 ft 
4u 


QQ 
— OO 


u 


4iy 


n fi 
— u.o, . 


. . , / u.o 




A 11 

All 


n 
U 


22 




1 


AO 

Uo 


on 
zU 





126 


—4.5 


72 


1 ftQQ I QQO 


Ground 


1 
1 


37 




3 


A A 
44 








893 


—0.1 


51.2 95 




A 11 

All 





zo 


x 


n 
u 


66 


—27 


n 

u 


^Q7 
oy f 


n 7^ 

— U. ( 0, 


. . . , yo 


1 ftQ7 I 
lOO / +0 ( 4 


Ground 


1 
1 


cc 



x 


2 


07 


— ZO 


u 


ouo 


n 9^ 

u.zo, 


fi4 
. . . , o^± 




A 11 

All 


n 
U 


26 


x 


u 


C^O 
OZ 


A 
4 


n 

u 


226 


— J..O, . 


. . , yo 


1 ft/l 1 1 Qon 

io4i-t-oyy 


Ground 


n 
U 


62 




3 


10 


Q 

— 


I 


877 


—0.2 


51.2 




A 11 
All 


U 


24 


x 


n 
u 


ft 1 
04 


1 A 
— 14 


n 

u 


QAn 

y^u 


2. 


fi4 


1 ft/i o_Lfion 
lo4z-t-0yu 


Ground 


9 
Z 


"1 7 


x 


2 


QQ 
OO 


Q 




n 

u 


ozo 


n 9^ 
— u.zo, 


(\A 
. . . , 04 




A 11 

All 


n 
U 


22 


x 


n 
u 


A 1 
4 I 


1 n 
— 1U 


u 


^79 
/ z 


— 1 25 


&n 

. . . , OU 


lOOz-hOyo 


Ground 


1 
1 


40 


x 


2 


no 

yo 


A 9 
— 4z 





509 


—0 4 

U.4, . 


51.2 




A 11 

All 


n 

u 


9^ 
zo 


x 


n 
u 


ftO 
Oz 


Q 1 
Ol 


u 


^9^ 
ozo 


1 

— 1, 


(\A 


I 7QO_l_ ^99 

I I oy-r-ozz 


Ground 


1 
1 


R9 
oz 


x 


Q 
O 


1 1 
11 


O 

— Z 


]^ 


814 


n n/L 

— U.U4, 


81 Q9 

. . . , 01 . yz 




A 11 

All 


U 


22 




f| 
U 


CO 
OO 


A 

— 4 


n 

u 


y 10 


n ^ 

— u.o, . 


fi4 


1 7/1 O I 7H1 

1 ( 4y+ / ui 


Ground 


1 
1 


83 


X 


2 


no 

yo 


70 
( z 





312 


-0.05, 


. .., 51.2 




A 11 

All 


n 
U 


25 


X 





77 

1 1 


7c; 
/ 





169 


-1.75, 


. . . , 56 


loUo-j- 1 o4 


Ground 


1 
1 


61 


X 


2 


yo 


90 

— zy 


1 


539 


-0.125 


. . . , 64 




All 





22 


X 





58 


-82 





928 


-0.8, . 


. . , 51.2 


1807+698 


Ground 


1 


61 


X 


2 


87 


-50 





567 


-0.07, 


. .., 71.68 




All 





28 


X 





46 


66 





320 


-0.1, . 


. . , 51.2 


1823+568 


Ground 


1 


66 


X 


2 


51 


-6 





817 


-0.3, . 


. . , 76.8 




All 





25 


X 





68 


16 





536 


-0.75, 


. . . , 96 


1828+487 


Ground 





95 


X 


1 


67 


-40 





631 


-0.28, 


. . . , 71.68 




All 





27 


X 





62 


-9 





379 


-0.9, . 


. . , 57.6 


1928+738 


Ground 





76 


X 


1 


59 


23 


1 


311 


-0.3, . 


. . , 76.8 




All 





40 


X 





54 


-42 





918 


-0.4, . 


. . , 51.2 


1954+513 


Ground 


1 


30 


X 


2 


54 


-3 





600 


-0.2, . 


. . , 51.2 




All 





25 


X 





48 


-34 





228 


-4, . . . 


, 64 


2021+614 


Ground 





99 


X 


1 


47 


66 


1 


112 


-0.4, . 


. . , 51.2 




All 





23 


X 





44 


-51 





401 


-1.2, . 


. . , 76.8 


2200+420 


Ground 


1 


57 


X 


3 


05 


-24 


2 


265 


-0.55, 


. . . , 70.4 




All 





22 


X 





46 


-27 


1 


239 


-1, ... 


,64 



Note. — Columns are as follows: (1) Source name; (2) VLBI baselines used 
in image; (3) FWHM dimensions of Gaussian restoring beam, in milliarcscc- 
onds; (4) Position angle of restoring beam, in degrees; (5) Peak intensity of 
image [Jy beam -1 ]; (6) Minimum and maximum contour levels, expressed as 
a percentage of peak intensity (intermediate contours are separated by factors 
of two). 
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Table 3 



Source Properties 



T ATT 
IAU 


vj trier 


vjpi. 




JVLBI 


i~>tot 


Name 


Name 


Trl 


z 


\ Tirl 


[ J yJ 


( L ) 


\ z ) 




(A\ 
I 4 ] 


\p) 


(b) 


0016+731 




Q 


1.781 


0.581 


0.8 


0133+476 


OC 457 


Q 


0.859 


1.985 


2.0 


0153+744 




Q 


2.338 


1.090 


1.1 


0212+735 




Q 


2.367 


3.045 


3.0 


0316+413 


3C 84 


RG 


0.017 


17.494 


22.2 


0711+356 


01 318 


Q 


1.620 


1.030 


1.0 


0814+425 


OJ 425 


BL 


0.245 


0.836 


1.0 


0836+710 


4C 71.07 


Q 


2.180 


1.762 


2.2 


0859+470 


4C 47.29 


Q 


1.462 


0.974 


1.3 


0906+430 


3C 216 


Q 


0.670 


0.684 


1.6 


0923+392 


4C 39.25 


Q 


0.699 


10.968 


11.0 


1624+416 


4C 41.32 


Q 


2.550 


0.668 


0.9 


1633+382 


4C 38.41 


Q 


1.807 


1.889 


2.4 


1637+574 


OS 562 


Q 


0.749 


0.703 


0.9 


1641+399 


3C 345 


Q 


0.595 


5.971 


8.2 


1642+690 


4C 69.21 


Q 


0.751 


0.878 


1.2 


1652+398 


MK 501 


BL 


0.033 


0.977 


1.7 


1739+522 


4C 51.37 


Q 


1.381 


1.933 


2.2 


1749+701 




BL 


0.770 


0.504 


0.6 


1803+784 




BL 


0.680 


2.040 


2.2 


1807+698 


3C 371 


BL 


0.050 


0.880 


1.6 


1823+568 


4C 56.27 


BL 


0.663 


1.047 


1.6 


1828+487 


3C 380 


Q 


0.692 


1.957 


5.3 


1928+738 


4C 73.18 


Q 


0.302 


3.324 


3.7 


1954+513 


OV 591 


Q 


1.223 


0.848 


1.4 


2021+614 


OW 637 


RG 


0.228 


2.664 


2.8 


2200+420 


BL Lac 


BL 


0.069 


4.069 


4.4 



Note. — Columns are as follows: (1) IAU Name; (2) 
Alternate name; (3) Optical identification, where Q = 
quasar, BL = BL Lacertae object, RG = radio galaxy; 
(4) Redshift; (5) Total 5 GHz cleaned flux density in 
VLBI image in Janskys; (6) Single dish 5 GHz flux den- 
sity at VLBI observation epoch in Janskys, estimated 
from the University of Michigan light curve. 
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Table 4 
Core Component Properties 



Source 


S 


Maj. 


Axial 


PA 


Fit 




Name 




axis 


ratio 


fdeg.l 

L o J 


Type 


r/io 12 ki 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0016+731 


0.118 


0.53 


0.45 


-64 


ZAR 




0133+476 


1.554 


0.17 


0.70 


-6 


ZAR 


3-9±i n o 


0153+744 


0.352 


1.17 


0.40 


-86 


Res 


0.034+°;°° 2 


0212+735 


1.206 


0.54 


0.72 


71 


Res 


0.28+Hs 


0316+413 


1.556 


1.17 


0.87 


62 


Res 


0.066+°°i 

—0.0.2 


0711+356 


0.078 


0.46 


0.19 


-36 


ZAR 


O.IO+A 1 ^ 

^ — 0.07 


0814+425 


0.448 


0.09 


0.53 


-66 


u 


5.4+l" f 

— 4.5 


0836+710 


0.603 


0.28 


0.46 


49 


ZAR 


0.83+A n j' 


0859+470 


0.447 


0.85 


0.16 


-31 


Res 


0.19+2-25 

—0.05 


0906+430 


0.483 


0.43 


0.15 


-34 


Res 


0.86+2? 

—0.2 


0923+392 


0.267 


1.01 


0.34 


-21 


n/a a 


0.038 


1624+416 


0.152 


0.37 


0.29 


31 


U 


0.20 + £? 


1633+382 


0.443 


0.15 


0.71 


-43 


Res 


1.46±g-i 


1637+574 


0.241 


0.13 


0.24 


27 


U 


3.13+' n 3 f 

— U . O 


1641+399 


0.811 


0.40 


0.63 


-26 


Res 


0.42+g-°5 

— U.U1 


1642+690 


0.396 


0.19 


0.14 


-22 


ZAR 


3.87+1"'' 


1652+398 


0.450 


0.24 


0.87 


32 


Res 


0.47+.2;;? 


1739+522 


1.745 


0.33 


0.88 


-29 


Res 


0.96+g-gl 

— U.UO 


1749+701 


0.264 


0.32 


0.52 


-50 


ZAR 


0.26+A 1 ! 

— U.UO 


1803+784 


1.374 


0.29 


0.63 


88 


Res 


1.28+21 

—0.2 


1807+698 


0.387 


0.28 


0.29 


71 


ZAR 


0.90+n'i 

— 0.3 


1823+568 


0.791 


0.60 


0.20 


16 


Res 


0.54^;} 


1828+487 


0.510 


0.42 


0.12 


-32 


ZAR 


1 o+inf 


1928+738 


0.416 


0.35 


0.34 


-21 


ZAR 


0.49±^| 


1954+513 


0.313 


0.40 


0.23 


-38 


ZAR 




2021+614 


0.078 


0.70 


0.40 


60 


n/a a 


0.02 


2200+420 


1.552 


0.42 


0.44 


10 


Res 


0.99t°i 



Note. — Columns are as follows: (1) Source name; (2) Fitted 
5 GHz flux density in Janskys; (3) Major axis of fitted Gaussian 
(FWHM) in milliarcseconds; (4) Axial ratio of component; (5) 
Position angle of component's major axis; (6) Visibility data con- 
sistent within errors with: U = unresolved component; ZAR = 
component having zero-axial ratio; Res = component that is re- 
solved on longest space baselines; (7) Brightness temperature of 
best component fit in units of 10 12 K , with range of possible values 
that fit the visibility data. 

a Componcnt too weak to determine errors in model fit parame- 
ters. 
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Table 5 
Jet Component Properties 



Cpt. 


r 


9 


S 


Maj. 


Axial 


PA 


log T b 




Name 


[mas] 


[deg.] 


[Jy] 


Axis 


Ratio 


[deg.] 


[K] 


Ref. 


( L ) 




C3~l 


V 1 ) 






(7\ 




rqi 


0016+731 


















C2 


0.65 


130 


0.254 


0.72 


0.37 


83 


10.9 




CI 


0.89 


137 


0.183 


0.95 


0.35 


-64 


10.5 




0133+476 




















o 

U.OO 


ri 
-Ol 


n 991 


u.oy 


n k 
U. ID 


u 


11 .u 






1 44 




n 1 09 

U. J.UZ 


u.oo 


44 


-Ul 


10 9 




PI 

V_/l 


1 7^ 

o. / o 




o 078 

U.U 1 O 


o.yo 


u. / O 


AA 


O.D 




0153+744 


















TZ> 

r> 


in IK 
1U. ID 






i.yi 


U.Oo 


-Dz 


1 A O 

lU.z 


i 


0212+735 


















C6 


0.34 


124 


0.880 


0.32 


0.37 


-7 


12.1 




L/O 


n fin 
U.DU 


117 
11 1 


n Q/i 7 

U.o4 ( 


n /i9 

U.4Z 


n nn 
U.UU 


99 
-ZZ 






PA 


2 21 




fl 3A7 


1 .DO 


fl A^ 


ss 

on 


±U. 1 




C3 


4.55 


107 


0.047 


1.70 


0.29 


89 


9.4 




C2 


7.00 


103 


0.050 


1.55 


0.41 


-52 


9.4 




CI 


14.06 


92 


0.179 


2.40 


0.68 


-53 


9.4 




0316+413 


















P9 


1 07 

l.U / 


-104 


n ^01 

U. Jul 


n or-; 
u.yu 


U.OU 


21 


1 Q 

iu. y 




Cl 


1.30 


-164 


1.109 


0.84 


0.67 


43 


11.1 




0711+356 


















P9 


n on 
u.yu 


9n 
-zu 


n nfii 

U.UDl 


n fii 

U.Dl 


n 7Q 
u. / y 


A^ 


1 n n 

1U.U 




pi 


^ sn 
o.ou 


93 
-zo 


n 7&a 
U. I o*± 


n S3 

U.OO 


n 7fi 
u. / D 


^3 
-Do 


in n 
iu.y 




0814+425 


















C4 


0.36 


100 


0.064 


0.27 


0.22 


-31 


11.3 




r<3 


1 nc; 


S7 

O 1 
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Note. — Columns are as follows: (1) Component name; (2) Distance 
from core in milliarcseconds; (3) Position angle with respect to core; (4) Flux 
density in Jy; (5) Major axis of fitted component in milliarcseconds; (6) Axial 
ratio of fitted component; (7) Position angle of component's major axis; (8) 
log of Gaussian brightness temperature in K; (9) Reference for component 
identifications: (1) Hummel ct al. 1997; (2) Paragi et al. 2000; (3) Lister & 
Smith 2000; (4) Klare et al. 2000; (5) Kameno ct al. 2000; (6) Murphy ct al. 
2000; (7) Tschagcr et al. 2000. 
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Fig. 1. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted images using 
all baselines. Angular scales are in milliarcscconds. 



Pearson-Readhead Survey From Space. I. 



13 




50/h pc 
0212+735: Ground Baselines 




I I D 




1/h pc 'J\ 


^^nM C2 


C1 






8F : 






\\\\\ 

% 


\ — yyw * i - 
' / 0316+413 [3C 84]: 




Ground Baselines 


10 5 


-5 -10 



50/h pc 




0711+356: Ground Baselines 

-5 



1/h pc 




0316+413 [3C 84]: Ground+Space Baselines 



50/h pc 




C2 
D 




0711+356: Ground+Space Baselines 

5 0-5 



Fig. 2. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 3. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 4. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 5. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 6. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 7. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 8. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 




Fig. 9. — Left column: naturally-weighted images using ground-to-ground baselines only. Right column: uniformly-weighted 
images using all baselines. Angular scales are in milliarcseconds. 
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Fig. 10. — Upper envelope of visibility amplitude distribution plotted against projected baseline separation for strong (left panel) and 
weak (right panel) sources in Pearson-Readhead sample. The flattening on baselines longer than an Earth diameter is indicative of strong 
unresolved core components. 
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Fig. 11. — Fitted Gaussian component size plotted against projected separation from the core component, for jets in the Pearson-Readhead 
survey. 
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Fig. 12. — Plots showing typical (u,v) plane coverages for VSOP observations using HALCA and the VLBA (left panel) and the EVN 
(right panel). 
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FlG. 13. — Simulations using complex core-jet source model (cj-complex) and jet PA — 60°. The images have a pixel size of 0.08 milliarcscconds. 
Panels arc as follows: (a) original source model; (b) (tt, v) coverage; (c) simulated HALCA image; (d) SNR ratio of each pixel in the difference image 
versus SNR in simulated image; (e) difference image; (f) Percentage error of each pixel in simulated image plotted against SNR ratio. The contours 
in panels (a) and (c) are 0.5,1,2,4,8,16,32, and 64% of the peak for both positive and negative contours. The contours in panel (e) are 8,16,32, and 
64% of the peak for both the positive and negative contours. The pixels clustered along the straight lines in panels (d) and (f) are located mainly 
off-source and have 100% error, since the model is essentially zero at these points. 
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FlG. 14. — Simulations using simple core-jet source model (cj-simple) and jet PA — 60°. The images have a pixel size of 0.08 milliarcscconds. 
Panels arc as follows: (a) original source model; (b) (tt, v) coverage; (c) simulated HALCA image; (d) SNR ratio of each pixel in the difference image 
versus SNR in simulated image; (e) difference image; (f) Percentage error of each pixel in simulated image plotted against SNR ratio. The contours 
in panels (a) and (c) are 0.5,1,2,4,8,16,32, and 64% of the peak for both positive and negative contours. The contours in panel (e) are 8,16,32, and 
64% of the peak for both the positive and negative contours. The pixels clustered along the straight lines in panels (d) and (f) are located mainly 
off-source and have 100% error, since the model is essentially zero at these points. 
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Fig. 15. — Dynamic range as a function of source position angle (PA) for the complex core-jet source model (cj- complex). The upper curve 
(labeled -D.Rrms) represents the ratio of maximum flux density in the image to the off-source RMS noise level. The lower curve (labeled 
DiJdiff) represents the maximum in the image divided by the maximum value (whether positive or negative) in the difference image. The 
major axis of the (u, v) coverage is held constant at PA= 0° for all models. 
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Fig. 16. — Dynamic range as a function of source position angle (PA) for the simple core-jet source model (cj-simple). The upper curve 
(labeled -D-Rrms) represents the ratio of maximum flux density in the image to the off-source RMS noise level. The lower curve (labeled 
DRdiff) represents the maximum in the image divided by the maximum value (whether positive or negative) in the difference image. The 
major axis of the (u, v) coverage is held constant at PA= 0° for all models. 



